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Résumé 

Comportement asymptotique des polynômes orthogonaux associés à un 
poids ayant un zéro fractionnaire sur le cercle. Applications aux valeurs 
propres de certaines matrices aléatoires unitaires 

Cet article s'intéresse aux polynômes orthogonaux correspondant aux poids de type 
(1 — cos0)"c où c est une fonction suffisamment régulière et a e] — i, i[. On décrit le 
comportement asymptotique des coefficients et de la valeur en 1 de ces polynômes et de 
leur dérivées lorsque N tend vers l'infini. Ceci nous permet d'obtenir une asymptotique 
du noyau de Christofel-Darboux associé à un tel poids et de calculer la loi conjointe 
des valeurs propres d'une famille de matrices aléatoires unitaires. Les démonstrations 
des résultats relatifs aux polynômes orthogonaux sont essentiellement basées sur les 
propriétés de l'inverse de la matrice de Toeplitz de symbole /. 

Abstract 

Asymptotic behavior of orthogonal polynomials on the circle, with res- 
pect to a weight having a fractional zéro on the torus. Applications to the 
eigenvalues of certain unitary random matrices. 

This paper is devoted to the orthogonal polynomial on the circle, with respect to a 
weight of type / = (1— cos 9)°'c where c is a sufficiently smooth function and a €] — ^, ^ [. 

We obtain an asymptotic expansion of the coefficients of this polynomial and of «l'j^-' (1) 
for ail integer p. Thèse results allow us to obtain an asymptotic expansion of the as- 
sociated Christofel-Darboux kernel, and to compute the distribution of the eigenvalues 
of a family of random unitary matrices. The proof of the resuts related with the ortho- 
gonal polynomials are essentialy based on the inversion of Toeplitz matrice associated 
to the symbol /. 
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1 Introduction. 



Il est bien connu qu'une matrice aléatoire est caractérisée par la loi de probabilité 
de la distribution de ses valeurs propres. Lorsqu'il s'agit d'étudier des matrices aléatoires 
unitaires, le cas le plus étudié est celui de l'ensemble circulaire de Dyson oii cette densité 
est de la forme, pour une matrice unitaire N x N (voir [5], [1], [6], [23]) : 

Pn{0i,---On) = cn n -e''N'- 

Il faut remarquer que dans le cas des matrices aléatoires unitaires c'est bien sûr la distibu- 
tion des arguments des valeurs propres qu'on étudie. 

Cet ensemble est considéré comme une bonne représentation des phénomènes physiques qui 
ont motivé initialement l'étude des matrices aléatoires. Il est naturel de chercher à l'étendre 
à un ensemble plus vaste oii la loi des valeurs propres est de la forme ( voir [T3] et [22]) : 

PNi9i,---9N)= n /(^J-) n (1) 

l<i<JV l<j<k<N 

(ensemble circulaire de Dyson généralisé) oii / est une fonction positive, intégrable sur le 
tore. Alors que, dans le cas de l'ensemble circulaire de Dyson, la fonction de corrélation des 
arguments des valeurs propres se calcule en fonction des polynômes orthogonaux associés 
à la mesure de Haar du tore ( [2Q] ) , dans le cas de l'ensemble circulaire de Dyson généralisé 
cette fonction de corrélation se calcule au moyen des noyaux de Christoffel-Darboux Kjsf 
associés aux polynômes orthogonaux unitaires ($„) oîi 1 < n < )associés au poids /. 
On définit généralement K]\f par (voir [19j) : 

, , ^-1 1 

KJv(e*^e*'') = Jf{e^')Jf{e^n E l-^m{e^')^m{e''') 

avec 

je 



On sait que, pour z / 0, on peut écrire aussi, si 9 ^ 6', et avec ^*i^{z) = z^^n (^) , 



'y ' (i_ei(^'-e)) 

Dans le cas unitaire, le lien entre ce noyau et la fonction de corrélation est rappelé dans 
le cours de la démonstration. La fonction de corrélation étant définie, un des problèmes 
relatifs aux matrices aléatoires est de déterminer la distribution des valeurs propres lorsque 
l'entier N tend vers l'infini. Ce problème se subdivise en fait en deux. D'une part il s'agit de 
déterminer la distribution des arguments de ces valeurs propres sur un intervalle ]e, 27r — e[ 
avec e > pour les matrices aléatoires unitaires ou la distribution des valeurs propres elles- 
mêmes dans l'intervalle ]0, -|-oo[ pour les matrices aléatoires hermitiennes (voir [21], [22], 
[23]). D'autre part on recherche cette distribution au voisinage de zéro pour les matrices 
aléatoires unitaires, ou au voisinage de zéro ou de plus l'infini pour les matrices aléatoires 
hermitiennes. 
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Dans le cas de matrices aléatoires Hermitiennes Gaussiennes et grâce à des résultats qui 
ont été obtenus sur les polynômes orthogonaux qui leur sont associés (on peut consulter, 
par exemple, [T], [3], [TT], [15] . . . ); on dispose de nombreux renseignements sur la loi des 
valeurs propres " au bord" du spectre (c'est à dire les valeurs propres proches de zéro ou 
tendant vers l'infini). 

Dans ce cadre Kuijlaars et Vanlessen [TT] étudient des matrices aléatoires Gaussiennes 
(hermitiennes) dont la loi des valeurs propres est donnée par 



1 ^ 



j=l i<j 

OÙ W]\f un poids de la forme Wn{x) = \x\'^°'e~^^^^\ pourx G M et F une fonction poten- 
tielle qui peut par exemple être une fonction polynôme. Les auteurs de l'article montrent 
que la fonction de corrélation des valeurs propres autour de l'origine est donnée par un 
noyau de la forme 

^a, ^ ^ ^Ja+l/2{'^u)Ja-l/2{T^v) - Ja-l/2i-ïï-u)Ja+l/2{'^v) 
J^(U, V) = TT^/U^/V ^ , 

2[u — V) 

on les fonctions Ja+1/2 ^t Ja-1/2 sont les fonctions de Bessel usuelles. 
Dans le présent article, nous obtenons un résultat du même type relativement à une classe 
de matrices aléatoires unitaires. En fait nous étudions les coefficients des polynômes ortho- 
gonaux de degré associés au symbole 

/(e*^) = |l-e'^|2°c(e*^) (2) 

avec c une fonction régulière sur le tore et a G] — ^, En effet nous obtenons dans un 
premier temps une expression asymptotique quand N tens vers l'infini des coefficients des 
polynômes orthogonaux puis nous en déduisons la loi des arguments des valeurs propres au 
voisinage de l'origine quand tend vers l'infini. 

Nous donnons d'abord (théorèmes [H et E]) l'asymptotique quand N tend vers l'infini des 
coefficients des polynômes et <I>jv. Pour ce faire nous utilisons des méthodes basées 

N 

sur des résultats antérieurs (se rapporter à [16] ou [I?]). Si l'on pose $Ar(e*^) = /3m e™^ 

notre énoncé se subdivise en trois parties qui correspondent à trois zones distinctes de 
l'ensemble des indices. Nous donnons d'abord les indices qui correspondent "aux bords" de 

XL U 

l'ensemble des indices, c'est à dire des entiers u tels que lim — - = ou lim — = 1 

et les indices qui se trouvent "au coeur" de l'ensemble des indices, ce qui correspond aux 

11 

entiers u tels que lim ~ avec < x < 1. Ces expressions nous permettent ensuite de 

calculer l'asymptotique de $^■'(1) et ^>^^^(1) (voir les théorèmes [î] et S]) , pour un entier p 
fixé. Nous pouvons alors calculer une représentation asymptotique du noyau de Christoffel- 
Darboux d'ordre N associé au symbole 9 \1 — e*^p"c(e*^). Si — vr < u,v < tt et u v, on 
a 

f iu/N iv/N^ N{uvY p{a,-u)p{a,v) - p{a,u)p{a,-v)e<''-''^ 

Aiv(e ,e ' ) = 0/ N — 7TT 7 N \-o[I\) [à) 

r"=(a)ci(l) i[u — V) 
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où p est une fonction qui sera précisée plus loin (voir les théorèmes [5] et [6|) , et dont l'ex- 
pression est différente selon que a G] — ^, 0[ ou que a g]0, \ [. 

Notons maintenant 9i , ■ ■ ■ 0^ les arguments contenus dans l'intervalle ] — vr, vr] d'une ma- 
trice aléatoire Un de taille N x N suivant la loi donnée dans l'égalité ^ avec un poids 

N 

W comme précisé en ([2]). Si l'on pose F^r = > ôgN les résultats précédents permettent 

i=l 

d'établir, si — vr < u < v < n (voir les théorèmes [5] et [6]) 
i) 

oià Fao est une fonction définie au moyen du déterminant d'un opérateur de Fredholm 
de noyau K = limAr^+oo -^f-- 
ii) Si q est un entier supérieur ou égal à 2 

(Xi V \ 

Heuristiquement, cela signifie que lorsque N tend vers l'infini il existe des valeurs propres de 
la forme e*^^ avec 9^ = 0{j^), mais qu'il n'en existe pas de la forme e^^^ avec = 0{^) 
pour q > 1. 

La première partie de nos résultats se rapproche en le complétant dutravail de Martinez- 
Finkelstein, Mac-Laughin et Saff qui dans [13] étudient une famille de polynômes orthogo- 
naux <î>7v orthogonaux pour le poids 

m 

W{z) = w{z)Yl\z-ak\'^^^ zeT, (4) 
fc=i 

où la fonction w est suffisamment régulière et où W vérifie les conditions de Szegô 

/ W{z)dz > et / ln{W{z))dz > -oo. (5) 

D ans l'article ([1^) ces auteurs donnent notamment le comportement asymptotique de 
$Ar(z) pour z ak 1 < k < m, et étudient également les zéros des polynômes $Ar, leur 
principal outil de démonstration étant la méthode de Riemann-Hilbert. 
A la fin de cet article, nous proposons un appendice dans lequel nous donnons des tableaux 
de valeurs numériques permettant de comparer <î>Ar(l) et les quantités données dans les 
théorèmes El et m pour des valeurs de a égales à —0,275, —0,150, 0,025, 0,1, 0,225. 



2 Notations et définitions. Enoncé des résultats 

Avant d'énoncer nos principaux résultats, nous introduisons des notations et définitions. 



2.1 Notations. Définitions 

Etant donnée / G L^(T) une fonction définie sur le cercle unité, on appelle matrice de 
Toeplitz de taille (A^ + 1) x (A^ + 1) associée au symbole / la matrice 

TN{f) = (f{i-j)^^..^^^ (6) 

V J 0<i,]<N 
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ou 

désigne la suite des coefficients de Fourier de /. 

Le lien entre les matrices de Toeplitz et les polynômes orthogonaux est bien connu (voir, 
par exemple, [12]). Ainsi nous avons 



Les polynômes y J sont souvent appelés polynômes prédicteurs. Comme le 

montre la relation ci dessus, leurs coefficients sont, à normalisation près, ceux de la première 
colonne de la matrice TN{f)~^ ■ Ce sont d'abord ces coefficients que nous donnons dans les 
théorèmes [T] et [2J 
D'autre part la relation 

mz) = z^'^Ni-), (9) 

z 

implique que les coefficients des polynômes orthogonaux correspondent, à une normalisation 
près, aux termes de la dernière colonne de la matrice (/). Nous utiliserons la notation, 
pour tout k dans Z 

/ : e'' ^ xHO) = e'"'. 

On notera (resp.iî_) le sous-espace fermé de L^(T) constitué des fonctions u dont les 
coefficients de Fourier f{k) = quand k < (resp. A; > 0). On notera par 7r± : L^(T) ^ H± 
les projecteurs orthogonaux respectifs. 

Nous rappelons la définition de l'espace de Beurling. On commence par introduire la notion 
de poids de Beurling : une suite {fij)jçz est dite poids de Beurling si elle satisfait aux trois 
propriétés suivantes 

' (i) fij > 1, j G Z 

(a) fij = fi^j, j e z (10) 

(m) fij+k < IJ'jlJ'k. j, k ez. 
La classe de Beurling associée au poids /x = (fij) est alors définie comme suit 

W^ = {we L\T) : wij) \< oo}. (11) 

Soit / un symbole tel que 

/(e^^) =1 1 - e*^ p'^ c(e*^) 

oii c est une fonction supposée être strictement positive sur cercle unité et appartenant à 
un espace de Beurling. Nous noterons Qa G la fonction analytique dans le disque unité 
ouvert, dite fonction extérieure, vérifiant 

/ =1 9a 

et Cl la fonction de vérifiant 

c = Cl Cl 

nous noterons par (/9^"^) la suite des coefficients de Fourier de l/ga, et enfin pour alléger 
les notations nous supposerons {(3^^) = 1. 
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2.2 Enoncé des résultats. 

2.2.1 Polynômes orthogonaux 

Nous nous intéressons, en premier lieu, au comportement asymptotique de la première 
colonne de l'inverse d'une matrice de Toeplitz Ti\i{f). Nous donnons d' abord un énoncé 
qui porte sur les coefficients situés sur les " bords " de la première colonne de l'inverse (voir 

m). 



Théorème 1 Soit Tj\f{f) une matrice de Toeplitz dont le symbole f s'écrit 

/=l 1-Xl'"c 

où a G M, \ a \ < ^ et où la fonction c > est supposée appartenir à une classe de Beurling 
avec jjLj > I pour tout entier j. Pour tout entier naturel k, tel que —>■ lorsque 
N ^ oo on a alors 
^) 

TNif),l,,, = - N^oo (12) 

avec Rn = 0{^^^) si < a < ^, et avec Rn = 0{-^) si — ^ < a < 0. 
ii) 

to'(/))„«_,,=(|§/3r")^+o(^). iv^». (13) 

La convergence est uniforme sur tout intervalle [0, [Ne]] et sur {N,N — [Ne]) pour e assez 
petit. 

On a d'autre part l'énoncé suivant qui porte sur les éléments du "coeur" de la première 
colonne de l'inverse ( |16j et [17j ) 

Théorème 2 Si Q < x < 1 est un réel donné, on a, avec les mêmes hypothèses que pour 
le théorèmeUl 

ci(l)(r^i(/))[^^,j^^_^ = i^„(x)iV--i + o(iV"-i), N^oo (14) 

où 

K^{x) = -^x'^~\l-xr. 
r(a) 

Le théorème\^la convergence est uniforme sur tout intervalle [ôi,Ô2] pour < ôi < §2 < 1- 

Avec les hypothèses que nous nous sommes données, les termes calculés ci-dessus corres- 
pondent aux coefficients de Ceux de en sont alors déduits à partir de la relation 
[9l Ce sont ces coefficients que nous donnons maintenant. 

N+l 

Corollaire 1 Posons ^]\f{z) = PuZ^ le polynôme orthogonal associé au symbole f = 

u=Q 

|1 — xP"c. Avec les mêmes hypothèses que pour les théorèmes{^ et{^ on a, si k/N 
quand N +oo 
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avec Rn = 0{^^^^) si < a < ^, et avec Rn = 0{^) si — ^ < a < 0. 



2. 



De plus, si < X < 1 est un réel donné, on a alors 

P[Nx] = Ka{l - x)N^-^ + o(iV"-i), iV ^ CX) (15) 

Comme nous l'avons rappelé dans l'introduction, le calcul du noyau de Christoffel-Darboux 
que nous nous proposons d'effectuer pour obtenir les noyaux annoncés nous oblige à déterminer 
avec précision les quantités $^^(1) et (<î>|^)^^^ (1) pour tout entier naturel p. C'est ce que 
nous faisons dans les théorèmes [3] et IH II est à noter que ces résultats, bien que donnés 
seulement dans le cas d'une seule singularité, prolongent les résultats de jl3J, puisque dans 
cet article les valeurs prises en les singularités par les polynômes orthogonaux sont inacces- 
sibles. 

Théorème 3 Avec les hypothèses des théorèmes ([ip et ^ on obtient 
-Si^<a<Oona 

= t;^^ / X--' ((1 -xr-i)dx+-)+ o(iv"). (16) 

r(a) ci(l) \Jo aj 
- D 'autre part, pour un entier j strictement positif fixé, dans le cas ^ < a < et pour 



un entier j positif fixé dans le cas < a < l 



on a 



($^)(J^) (1) = / x'^-^+^il - x)^dx + o{N^+^). (17) 



r(a) ci(l) 7o 

Remarque 1 On peut remarquer que pour un entier j strictement positif et suffisamment 
petit devant N toujours avec l'hypothèse — ^<a<^ on a en fait 

($^)(^-) (1) = jV^+.-r(a + j)r(a + l) 1 

^ ^' ^ ' r(2a + i + l) r(a)ci(l) ^ ^ 

On obtient le même genre d'énoncé pour les polynômes orthogonaux. 

Théorème 4 Avec les hypothèses des théorèmesUl etlM on obtient 

- Si < a < et si j un entier positif on a 

= CfïTÏÎ) - - 1) + ^) + °(^"")- (18) 

- D 'autre part, pour un entier j positif fixé dans le cas < a < ^ nous pouvons écrire 

{^N^^^ (1) = -pTT-m / ^""'(1 - ^)"^'^^ + o{N^^'). (19) 
T[a) ci(l) Jq 
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Remarque 2 On peut là aussi remarquer que pour un entier j strictement positif suffi- 
samment petit devant N et toujours avec ^ > a > on a 

^ ^ > r(2a + i + l) r(a)ci(l)^ ^ ^ 

2.2.2 Matrices aléatoires unitaires. 

Soit maintenant Un une matrice aléatoire unitaire de taille A^. Notons 9^ 1 < i < A'^ les 
arguments des valeurs propres de Un que l'on suppose classés dans l'ordre croissant dans 
] — TT, +7r]. On suppose que ces arguments ont pour densité de probabilité 

N N 

avec /(e*^) = |1 — e*^p"ci(e*^) et oii ci est dans une classe de Beurling d'indice supérieur 

N 

OU égal à |. Nous posons ensuite Fn = ''^^^qn où Ôqn désigne la mesure de Dirac de 

support . Dans la suite nous supposerons enfin que u et v sont deux réels tels que 
— TT < u < V << TT. On a alors les énoncés suivants qui précisent la loi au voisinage de 1 
des valeurs propres des matrices aléatoires unitaires étudiées dans cet article. 

Théorème 5 Soit < a < ^. Posons, pour tout n G M 

Nous pouvons alors écrire, avec les hypothèses ci-dessus, et pour tout entier positif m 

OÙ I = [u,v]et où K, est l'opérateur de Fredhom de noyau K avec 

- si 9^9' 

(f) f)>) f -0)^{a, 9') - ïp{a, 9)i){a, -9')e'^' \ 

^ r2(a)c(l)e-*^ V i{9' - 9) )' 

- et si 9 = 9', 

Q2a 



(-25JÎ (V'(a, 9)T{a, -9)) + \ip{a, 9) 



2^ 



T\a)c{l) ■ 

Remarque 3 On a en particulier le résultat 

^r5îoo^(^^^l^'^] = °) ^^^^ 
Le théorème suivant permet d'envisager le cas oii a. est négatif. 
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Théorème 6 Soit — i < a < 0. Posons, pour tout u,v ^ 



1 

'iP{a,u)= I x"-^((l-x)"e-*"^-l)dx + -. 

Jo " 



Nous pouvons alors écrire, avec les hypothèses du théorème\^ 

lim p(F^[^-]=m) = ^-J^{-) det[(Id-7/C)|i.(,)] (22) 
N~^+oo V iV iV / m! \a7/ ^ ' 7=1 



où I = [u,v] où te est l'opérateur de Fredhom de noyau K avec 

- sie^e' 

nf^ / e'^tp{a, -e)tp{a, d') - tP{a, 9)tp{a, -d')e'' 

[0,0) 



25R (î>{a, 0)T{a, -9) ) + \î}{a, e)\ 



r2(a)c(l)e-^^ V i{e' - 9) 

Ou encore si 9 = 9' , 

r2(a)c(l)' 

Remarque 4 On a en toujours le cas particulier intéressant suivant 

Fn[j;^, ^] = 0) = det[(/d - /C)|i2(,)] (23) 

Théorème 7 Pour tous entiers naturels q et m vérifiant q > 1 et N > m > 1 et pour tout 
réel a tel que ^ < a < ^ on a, toujours avec les mêmes hypothèses qu 'au théorème 

lim p(FN\—,—]=m)=0(N^^-P'>). (24) 

3 Preuves des résultats 

3.1 Preuve des théorèmes [T] et [2] 

Les points i) et ii) du théorème [T] ont été démontrés dans La preuve est basée sur 
une formule explicite de l'inverse que nous avons obtenue dans un précédent travail ( [IGJ 
ou [IZ]). Nous pouvons en effet écrire pour tout entier A:, < A; < iV + 1, l'expression 



^ (E (^u^^N,a {{k - U)/N)^ ^ (1 + 0(1)) (25) 



avec 

m=0 



TT 



oii on a posé, pour z G [0, 1] 

_p(°) = V = V = + V = + 

2m,N\ ) jo + ^Jl+k + N^'" . ^ n hm-l + j2m-2 + (26) 

JO J1=0 J2m-1=0 ^ ^ 

1 Y^OO 1 1 

^i2m=0 j2rn+j2 m— 1 
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où Na = N + 1 + a. Dans ([9j), Inoue montre que ;^-^2m iv('-') approximé par une 
intégrale dont la valeur est a^. En combinant ce résultat à un théorème d'approximation 
de Bleher (pj) permettant d'estimer l'écart entre l'intégrale et la somme, on obtient le 
développement 



N ^iV2- 



m=0 

Ces résultats permettent alors d'obtenir le point i) du théorème. Le point ii) est une 
conséquence de la relation de récurrence de Szegô (voir [7]) : pour | z |= 1 



OÙ jN = — ^n+i(0) est appelé coefficient de Verblunsky. L'uniformité est une conséquence 
de l'uniformité de la formule [25l 

Pour le théorème El le lecteur se reportera à [16j ou à [17]. 

Dans le calcul précédent nous avons déterminé le comportement asymptotique du polynôme 

prédicteur ^*j^{x) = ^ (■^))fc+i,i ^ déduit alors celui de $jv via l'opération 
k=o yTNif))î^ 



z 



Au lieu de considérer <î>jv(e*^) nous considérons ^*^{e^^) le polynôme prédicteur as- 



N 



3.2 Preuve du théorème [3] 

socié. La relation [9] montre que $Ar(l) = ^^(1). Posons ^^(e*^) = ^ /37v,i(e*^)^ Ecrivons 

z=o 

[Ne] N-[Ne] N 

= Si + S2 + Ss avec 5*1 = '^i3n,i, 5*2 = ^ /3n,i, 83= ^ PN,h avec e 

(=0 l=[Ne]+l l=N-[Ne]+l 

un réel positif tendant vers zéro.. En utilisant le théorème [T] nous obtenons 

^(a+l) l 



^1 = E(^r'-ij^r''+oi 



(^+1) _ {a+2) _ W^rX — 



Comme on a d'autre part 



?{"+!) _ 

r(a + l)ci(l)" 



A = T^T— ^rTTTTT^" + OiH et 



"r(a + 2)ci(l)' 
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nous pouvons conclure 

r(a + l)ci(l) V a + 1 

ou encore 

^l = w ^^.^ [iV6]" + 0(iV"). 

r(a + l)ci(l) 

Considérons maintenant la somme 82- Cette fois ci d'après le théorème [2] on a 

5, V — (1-—)°. 

^ r(a)ci(l) iV ^ A^^ 



Nous allons maintenant distinguer pour notre démonstration les cas a > et a < 0. 
Le cas ^ < a < 0. 
Nous pouvons écrire 



1 



r(a)ci(l) 
1 



AT" I / ^ ((1 - t)" - 1) (it 

([Ne]+1)/N 



J([iVe]+l)/iV / 



— -3— -AT" ( ((1 - t)" -l)dt + - 

r(a)ci(l) VJo " 



ar(a)ci(l) 



Compte tenu de = O ( 1 pour A; G [A^ — [A^e]] on obtient évidemment comme plus 

haut 53 = 0(A^°e"+i) = o{N'^). D'où 

$^(1) = 5i + 52 + 53 = 'p^-^-—N- ( f ((1 - tr -l)dt+-\+ o{N-). 

r(a)ci(l) \Jq a) 

Le cas < a < ^. 

Puisque a > on a tout de suite 5*1 = 0{N°') et d'autre part 

«^ = r{i(ï)''°(i''°"''-''"*)^ 

Puisque l'on a toujours ^3 = o{N") on peut écrire 
Calcul des dérivées 

Démontrons maintenant la formule dans le cas où p est un entier strictement positif. Avec 
les mêmes notations que précédemment nous avons 

AT 

(cî>^)(î')(e*^) = ^/(/ - 1) . . . (/ _p + i)/3^,(e^ey^ 
i=P 
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[Ne] 

Nous sommes ramenés à calculer (^>|^)(^)(1) = Si + S2 + avec 5*1 = 'Y^PPn,^ 5*2 = 

1=0 

N-[Ne] N 

Ff3j\r^i, S3 = I^Pn,u avec e un réel positif tendant vers zéro.. En utilisant 

l=[Ne]+l l=N-[Ne]+l 

le théorème ([T]) nous obtenons 

([Ne] 2 , 

= 0(iV°+Pe"+P) = o(iV"+P). 

On obtient de même = o{N"''^p), et en utilisant le théorème ([2]) il est facile de vérifier 
que 

^2 = — J— -AT^+P ( [ - t)''dt] + o(Ar"+P), 

r(a)ci(i) V7o y 

ce qui achève la démonstration de ce cas. 



3.3 Démonstration du théorème [4] 



Remarquons tout d'abord que bien stir <I'Ar(l) = $^(1). Pour le calcul des dérivées 
d'ordre supérieur ou égal à un nous allons là aussi devoir distinguer le cas a positif du cas 
a négatif. Là aussi trois cas sont à distinguer. 



TV 



Supposons d'abord a positif. Considérons les sommes (3n-uu{u — l)(n — 2) • • • (n — p + 1) 

u=p 

N 

dont la partie principale est en fait P^-^u^- Là aussi nous pouvons utiliser une décomposition 

u=p 

N [Ne] N-[Ne] N 

Pn-uUP = f^N-uVP + Y Pn-uU^ + Y ^N-uUP, 
u=0 u=0 u=[Ne] N-[Ne] 

OÙ e est un réel strictement positif qui tend vers zéro. En utilisant le théorème [T] nous 
pouvons écrire 

[Ne] _ [TV.] 

et l'on vérifie facilement que ce dernier terme est d'ordre 0(A^"+Pe°'''"^^^) = o{N'^^p). Nous 
avons de même 

N [Ne] , 2- 



iv L^'cj , 2 \ 



TV-[TVe] u=0 
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Il est là encore clair que ce terme est d'ordre 0(A^"^^)e""^^) = o(A^"^^). Reste à calculer 

N-[Ne] N-[N€] 

Pn-uuP ce qui peut aussi s'écrire jïu{N — uY , ou encore avec le théorème [T] 

u=[Nt\ u=[Ne] 
N-[Ne] N-[Nt\ 

u=[Ne] y J i-\ J „=[jve] 

La formule d'Euler et Mac-Laurin et l'hypothèse e tend vers zéro permettent de conclure 

P4N-ur= / (i-tr+H'^-'dt+oiN'^+n, 

nfi^.] r(a)ci(l)7o 

et finalement 

Ce qui est l'énoncé attendu. 

Nous allons maintenant traiter le cas a négatif. 

En introduisant la même décomposition que dans le cas précédent, considérons tout d'abord 

N 

la somme Si = Pn-uU^- On a 

N-Ne 

N . 2 \ i^''] 



u=N-[Ne] ^ ^ «=0 ^ ^ 

fc=0 \«=0 ^ ^ 

^ ^ V (a + ^)r(a) (a + A: + l)r(a + 1) 7 ci (1) 

T^p+a a Tvrp+a a 

+ eO(iVP+"e") = — + o(iVP+"e"). 



r(a + l)ci(l) ' ' r(a + l)ci(l) 
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N~[N€] 

On a d'autre part , si l'on pose 5*2 = Pn-uI- 

u=[Ne] 

N-[N£] N-[Ne] 
u=[Ne\ u=[Ne\ 





fQ-l 


T{a] 


ICl(l) 


NP- 


fo-1 


T{a] 


ICl(l) 


NP- 


fo-1 


T{a: 


ICi(l) 


NP- 


fa-1 


T{a] 


ICl(l) 



0-1, 



(1+0(1)) 



N 

u=[Ne] 



fl-e 

-e T\Tp+a a AfP+a 
- t)'^+P - l)dt - ^^TT + 



r(a + l)ci(l) r(a + l)ci(l) 

La troisième somme intervenant dans la décomposition étant clairement négligeable, on a 
finalement 

^Ni^) = rlT'n, ( r^t^'Hil - tr^" - m + -)+ oiNP+^-'), 
r(a)ci(i) VA "/ 

ce qui donne immédiatement le résultat. 
3.4 Preuve du théorème [5l 

3.4.1 Déterminant des opérateurs de Fredholm dans L^(T) 
Pour {0, 9') dans ] — vr, vr] x] — vr, vr] nous posons 

N 

k{e,e') =Y,^rn{0)fn.{0') 
m=l 

011 pour tout entier m tles fonctions (pm et fm sont dans L^(T). Si X est un intervalle ] — vr, vr] 
on peut alors définir un opérateur K sur Lp'iX) en posant, pour tout x dans -L^(X), 

K^{e) = [ k{e,e')x{e')de'. 

Jx 

Il est alors connu (voir, par exemple, [8] ou [18] ) que det(/ + K) est la quantité définie par 

N ^ 

det(/ + K) = 1 + V — - / det{k{ei,ej)ij=i,... ^mdOi ■ ■ ■ dam 

OU encore en posant que l'intégrale sur vaut 1 on peut écrire, de manière plus synthétique 

N ^ 

det(/ + K)=y" — / det{k{ei,ej)i.j=i,.. ^mdOi ■ ■ ■ dOm 
^0 
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De même on peut définir la trace de K par 



TtK 



X 



k{e,e)de. 



D'autre part si (i^n)neN est une suite d'opérateurs ainsi définis qui converge vers un 
opérateur A au sens des opérateurs sur Lp'{X), on peut poser 

à(it{Id + A)= lim det(/ + i^„) et TtA= lim Tr(i^„). 

n— ►oo ri— >oo 

On sait que A est alors un opérateur fixé (voir [8j) la fonction A — > det(/ + A^) est une 
fonction entière sur C dont le développement en série entière est donné par 

det(Id + A^) = l + f;^^A". 



71=1 



Les coefficients c„ étant définis par 



Tt{A^) 



Cn{A) = det 



n — 1 
Tr(A) 





n-2 








\ 





Tr(^'^-i) 
V Tr(^") 

Ce qui peut s'écrire aussi 



Tr(A"-2) Tr(^"-3) 



Tr(^) 1 
Tt{A^) Tï{A) ) 



det(/d + A^) = exp ^ Tr(A"^)A'^ 



\m=\ 



ml 



On peut alors énoncer le résultat : 

Propriété 1 L'opérateur Id + A est inversible si et seulement si det{Id + A) ^ 0. 



3.4.2 Polynômes orthogonaux et fonction de corrélation 

Considérons l'égalité 



N 



if,* ( ™ 



Nous pouvons écrire, pour un entier fixé 



j=o -^^ 



(A^o + 1)! 



En effet la démonstration du théorème [3] permet d'obtenir , si j > /cQ) et en remarquant 
que /(Z -!)(/- 2) <V : 



ci(l)r(a) Vio 



1 -x)^ 



■dx + o(iV") 
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uniformément par rapport à j. Ce qui donne, toujours en utilisant le théorème [3] 

Ceci peut encore s'écrire 

r(a)ci(l)$^(e^^/^) = iV" (^j^ e^^^x"-i(l - x^dx + Riko)^ . (29) 

avec 

Si maintenant nous nous fixons un réel e positif, il est clair que l'on peut choisir le réel fco 
tel que |iî(A;o)| < e. Nous pouvons donc finalement écrire 



r(a)ci(l)$^(e^^/^) = iV" (^J^ e^^^x"-i(l - x)"dx + o(l)^ . 



(30) 



On obtient de même que 

r(a)ci(l)$^(e*^/^) = iV° (^j^ e*^(i-^)a;"-i(l - x)''dx + o(l)^ . (31) 
Introduisons maintenant le noyau de Christofel-Darboux d'ordre A^, à savoir la quantité 

771=0 

OÙ les hra sont des constantes de normalisations définies par 

r f{e'')<èmie''^)<èn{e'')de = h^5m,n- 

J —n 

On sait que l'on peut écrire alors, si 6 ^ 9' (voir [19]) 

En posant dans la formule ci-dessus 6 = j^,6 = jj on obtient, en utilisant les formules 
(El et (El 



'K{u,v) + o{N 



avec si n 7^ î; 



(32) 



K{u,v) = / e-™^x"-i(l -x)"(ixy e^^^x"-i(l - x)"dx - 



1 

e" 
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Ou en posant 



1 







K{u, v) = Tpia, —u)ip{a, v) — ip{a, u)'tp{a, —v)e^^'" 
et en remplaçant y/J et (1 — eH^*/^-'"/^)) par un équivalent, il vient, toujours si u ^ v, 

r2(a)ci(l)K^(e»/^,e»/^) = J- N (uv)" + o{N) . 

Pour connaître la valeur u = v considérons la dérivée au point f = u de la fonction Gu 
définie par 

Gu{v) = é'^^l){a, -u)ip{a, v) - e'''"ip{a, u)ip{a, -v). 

La limite quand v ^ uà& (e*"'^^, e™/^) est alors j^N{u)'^" Ce qui donne, 

en posant r(a,u) = — j;)"e*"^(ix, 

r2(a)ci(l)i^^(e™/^,e*"/^) = -^iVn^" (-2^{i^{a,u)T{a,-u)) + \i^{a,u)\^) + o{N). 

On sait d'autre part (voir |12j ) qu'avec les notations du corollaire [T] 

hm= r f{e'')\<^N{e'Yde, 

J —n 

c'est à dire 

Nous avons d'autre part obtenu dans un précédent travail ( [I6],[T7]) que 

Finalement l'hypothèse Pq = 1 nous donne, avec u ^ v, 

r2(a)ci(l)i^^(e^"/^, e^^/^) = iV(n^)"^i^ + o{N), 

t[u — V) 

et si M = f 

r2(a)ci(l)i^^(e*"/^,e^"/^) = Nu^" {2^{i;{a,u)T{a,u)) + \'iP{a,u)f) + o{N). 
Si maintenant on désigne par R^^^ la fonction de corrélation définie par 

iîf Vi, ^2, ■■■ ,0k) = j-^—^ [ Pn{Oi, • • • , eN)d9k+i ■■■deN 

[N - k)\ Jy^^^yN-k 

avec 

N 

Pn{9i, ... ,eN) = ^\l- e^'^^f -ci(0,) n W'^ - e^'f • 
On sait qu'on a alors la relation fondamentale 
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Propriété 2 pour tout intervalle I contenu dans ] — 7r,7r[ et tout entier m G {0, • • • ,N} 
on a 



P 



{i <N, OiE 1} 



m 



(-l)'" A (-1)'^ 



E 



m '■ — ' (k — m)\ Jjk 

k=m 



Rr>{9)d''{9). (33) 



Des calculs classiques sur les polynômes orthogonaux donnent la relation 

Rf\0i ■■■ek) = det[KN{ei, ej)ij=i...k]. 

C'est à dire que 



({i<N, 0.^1} =m) = ^^Y.ir^^l 

\ J m {k — my.Jjk 

k=m 



det[KN{ei,ej\,,j=i...k]d^{9). 



(34) 



On obtient donc finalement, en posant maintenant / = [u,v\, — tt < u < t; < tt et 7;v 
[:^> la formule 

Pi{i<N, OiE In} =m)= ^ / det[i^iv(^i, %]i,,=i...fe]d*^(^). 



En posant le changement de variables Uj = N9i on obtient 

p[\{i<N, 6iGlN}\='' 



(35) 
(36) 



N 



i-ir ^ (-1)' 



m — ' (k — m)\ jjk 

k=m ^ ' 



det [iriv(ui/iV, UjlN\j=x...u\d^{u), 



ou encore 



Py{i<N, e^GlN} =mj=^Yl (k-m)\ J k ^jkj=^-k]dHu), (37) 

k=m 

avec,si u ^ v 



Hn{u, v) 



(uv) r+o(l), 



r2(a)ci(l) i{u — v) 



et si u = V 



Hn{u, u) = y (-2SR (V'(a, «)T(a, -n)) + |V'(a, u)!') + o(l). 

En revenant à la définition des déterminants des opérateurs de Fredholm nous pouvons 
écrire 

^ (-1)'' r 

Y^-rh det[iï;v(«i,%)i,,=i...fc]ci'^(^^),= det[(/-i/jv)|L2(/)] (38) 
it=o ■ 

Ce qui donne, au moyen d'un changement de variables. 



J2 / det[HNiui,Uj)ij=i...k]d''{u), = det[(/ - 7Hn)\l^i)] 

k=o -^i" 



(39) 
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ou encore 



N 

E 

k=m 



(-1)'7' 



k-,k—m 



{k — m)\ J jk 



■det[(/-7i/iv)|L2(7)]. (40) 



soit 



TV 

E 

k=m 



(k — m)\ Jjk 



detl 



HN{Ui,Uj)ij=i...k] d^iu) = \^j;^det[{I -jHN)\L\I)]j (41) 



Et par passage à la limite on obtient finalement 



lim P 

7V^+oo 



{i < NOi G /iv} = m] 



) i \ det[(/d-7/C)|i.(,)] 



m! y ^7 



7=1 



(42) 



Ceci termine la démonstration du théorème [5l Le deuxième point du théorème est alors la 
traduction du cas m = 0. 

La démonstration du théorème [6] repose bien sûr sur les mêmes idées que celle ci. 

3.5 Démonstration du théorème [6] 

Nous pouvons de nouveau écrire 



^?.(e*') = E 



:cï>^)œ(l). 



Comme dans le cas a positif le théorème (l3|) permet donne alors pour j > 1 : 



c(i)r(a) 70 



-dx + oiN" 



Ce qui donne finalement, en utilisant les mêmes arguments que dans la démonstration 
précédente, 



r(a)ci(l)^. 



N' 



li^$^(e^^/^) = Y. 



1 - xTdx 



+ / ((1 -x)" - l)(ix + - + o(l) 

Jo « 

Tout ceci permet de conclure que 

r(Q)ci(l)$^(e^^/^) = N''( C x"-i(e^^^ - 1)(1 - x)°dx+ 



(43) 



+ C x'^-^iil-xT -l)dx + - + o{l)\ 
Jo a J 



ou encore 



r(a)ci(l)$^(e^^/^) = N'' (^j^ (e^^^(l - x)" - l) + ^ + o(i; 



(44) 
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On obtient de même, en utilisant la relation$^(z) = z^^]\f{^) que 
r(a)ci(l)$7v(e*^/^) = iV"e^^ (^j\e-'^''{l - x)" - l)x"-^dx + 1 + o(l)^ . (45) 



C'est à dire que 



Désignons maintenant par K]\f le noyau de Christoffel-Darboux d'ordre N. On pose de 
nouveau 6 = ^,et9' = et en on utilise les formules (j43]) et (f45]) . On obtient alors,par les 
mêmes méthodes que précédemment, et pour u ^ v 

r2(a)|ci(l)p^^(e»/A^,eWA^) = ^Ar^"^/(eW^)^/(eW^) ^(;;' +/_(^^^^^ (46) 
avec 

K{u, v) = tp{a, —u)xp{a, v) — 'ijj{a, u)'>J'{a, —v)e'^^'"~'^\ 

et 

1 

tljia,u)= / x°-i((l-x)"e-™^-l)dx + -. 

Jo « 

En remplaçant ^/f et (1 — e*("/^-^/^)) par un équivalent, il vient si u ^ v 
r2(a)ci(l)^;v(e^"/^,e-/^) = J-Ar(-H°^^i^ + o(iV). 

La valeur pour u = v s'obtient par un simple passage à la limite. Le reste de la démonstration 
est alors identique à celle du théorème [SJ 

3.6 Démonstration du théorème [7] 

Nous pouvons de nouveau écrire pour p > 1 

^ (je/NP 1 \j 

j=o ^■ 
Là aussi le théorème [3] permet d'obtenir : 

_^($^)0)(1) = / ^ V cix + o(iV"). 

i! V at; V ; ivj(p-i)ci(l)r(a) io j! ^ ^ 

Ce qui donne finalement 



Ce qui peut aussi s'écrire 

r(Q)ci(l)$^(e^^/^'') = N° (^j^ e*''^/^''"'x"-^(l - x)"dx + o(l)j . (47) 
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On obtient de même que 

r(a)ci(l)$^(e^^/^) = TV" ^ie{i-x)/N^~\a-i^^ _ ^^a^^ ^ ^^^^^ (4g) 

C'est à dire que 

r2(a)ci(l)i^iv(e™/^',e™/^') = ^N^^ \[J{é^)\[J{ë^) (49) 

K{u/NP^^,v/NP-^) + o{N'^'^) 

^ Y — çi{u/NP-v/NP) ■ 

Ce qui permet de calculer les probabilités P (^{i < N, 9i £ Inp} ^ = ^} pour tous les 
entiers m,p > 1. En posant I = [u, v], —n < u < v < ir et Inp = [7^, j^p], on a 



N 



P 



{i <N, OiG Inp} 



m 



(-ir ^^ (-1) 



m ^ — ' (k — m)l j jk 

k=m ^ ' •'^N 



àei[KN{ei,ej]i^j=i...k]ê{e). 



En posant le changement de variables Ui = N^Oi on obtient 
p{\{i<N, O.^In} 



(50) 



m ^ — ' {k — m)\ Jk 

k=m ^ ' 



N'P'' det[KN{ui/NP, Uj/NP],^,=^...j,]d''{u), 



ou encore 



P(\{i<N, 9,elN} 

N 



m 



m ^-^ {k — m)\ Jjk 



det[HN{ui/NP'\uj/NP-\,=i...k]d''{u), 

Puisque pour tout entier naturel k G [1, A^] la quantité det[HN{ui/N^^^ ,Uj /NP^^)ij=i...k] 
tend vers zéro quand tend vers l'infini par continuité de la fonction Hn il est clair que 
le résultat annoncé dans l'énoncé est acquis. 



4 Appendice, résultats numériques 

Dans chaque tableau la première colonne donne 'î>Ar(l) la deuxième colonne donne AN~'^ 
où A est une constante induite par les résultats des théorèmes [3] et [H Les valeurs de N 
étudiés sont entre 400 et 640. Rappelons que par rapport au corps de l'article d = —a. 
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4.0.1 d=-0.2750 





n{1) 




^iv(l) 




$iv(l) 


AN--^ 


$ 


iv(l) 


AN-'^ 


11 


.7692 


11, 


.7688 


1 11.8282 


11.8279 


11.8811 


11, 


.8809 


11 


.9334 


11, 


.9333 


11 


.7746 


11, 


.7742 


11.8335 


11.8332 


11.8863 


11 


.8861 


11 


.9386 


11, 


.9385 


11 


.7800 


11, 


.7796 


11.8388 


11.8385 


11.8916 


11, 


.8914 


11 


.9437 


11, 


.9437 


11 


.7854 


11, 


.7850 


11.8441 


11.8438 


11.8968 


11, 


.8967 


11 


.9489 


11, 


.9489 


11 


.7907 


11, 


.7904 


11.8494 


11.8491 


11.9021 


11, 


.9019 


11 


.9541 


11, 


.9541 


11 


.7961 


11 


.7958 


11.8547 


11.8544 


11.9073 


11, 


.9071 


11 


.9593 


11 


.9592 


11 


.8015 


11, 


.8011 


11.8600 


11.8597 


11.9125 


11, 


.9124 


11 


.9644 


11, 


.9644 


11 


.8068 


11, 


.8065 


11.8653 


11.8650 


11.9177 


11, 


.9176 


11 


.9696 


11, 


.9696 


11 


.8122 


11, 


.8118 


11.8705 


11.8703 


11.9229 


11, 


.9228 


11 


.9748 


11, 


.9747 


11 


.8175 


11, 


.8172 


11.8758 


11.8756 


11.9282 


11, 


.9281 


11 


.9799 


11, 


.9799 



1.8228 11.8225 
4.0.2 d= -0.1500 



$iv(l) 


AN-"^ 


$iv(l) 


AN-'^ 


^iv(l) 


AN-"^ 


$iv(l) 


AN-'^ 


5.4235 


5.4234 


5.4449 


5.4448 


5.4542 


5.4541 


5.4711 


5.4723 


5.4249 


5.4247 


5.4463 


5.4462 


5.4528 


5.4554 


5.4723 


5.4736 


5.4262 


5.4261 


5.4476 


5.4475 


5.4555 


5.4567 


5.4736 


5.4749 


5.4276 


5.4274 


5.4489 


5.4475 


5.4568 


5.4580 


5.4749 


5.4762 


5.4289 


5.4288 


5.4383 


5.4382 


5.4581 


5.4593 


5.4633 


5.4646 


5.4303 


5.4301 


5.4396 


5.4395 


5.4594 


5.4606 


5.4646 


5.4659 


5.4316 


5.4315 


5.4410 


5.4408 


5.4607 


5.4620 


5.4659 


5.4672 


5.4329 


5.4328 


5.4423 


5.4422 


5.4620 


5.4633 


5.4672 


5.4685 


5.4343 


5.4342 


5.4436 


5.4435 


5.4633 


5.4646 


5.4685 


5.4697 


5.4356 


5.4355 


5.4449 


5.4448 


5.4646 


5.4659 


5.4698 


5.4710 


5.4370 


5.4368 


5.4463 


5.4462 


5.4659 


5.4672 


5.4711 


5.4723 


5.4383 


5.4382 


5.4476 


5.4475 


5.4672 


5.4685 


5.4723 


5.4736 


5.4396 


5.4395 


5.4489 


5.4475 


5.4685 


5.4697 


5.4736 


5.4749 


5.4410 


5.4408 


5.4502 


5.4475 


5.4698 


5.4710 


5.4749 


5.4762 


5.4423 


5.4422 


5.4515 


5.4501 










5.4436 


5.4435 


5.4528 


5.4515 
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4.0.3 d= -0.0250 

$jv(l) AN-'^ 

2.3768 2.3768 

2.3769 2.3769 

2.3770 2.3770 

2.3771 2.3771 

2.3772 2.3772 

2.3773 2.3773 

2.3774 2.3774 

2.3775 2.3775 

2.3776 2.3776 

2.3777 2.3777 

2.3778 2.3778 

4.0.4 d= 0.100 

$jv(l) AN-'^ 

0.9706 0.9707 

0.9704 0.9705 

0.9702 0.9703 

0.9701 0.9702 

0.9699 0.9700 

0.9698 0.9699 

0.9696 0.9697 

0.9695 0.9695 

0.9693 0.9694 

0.9691 0.9692 

0.9690 0.9691 

4.0.5 d= 0.2250 

$iv(l) AN-'^ 

0.3566 0.3568 

0.3565 0.3567 

0.3563 0.3565 

0.3562 0.3564 

0.3561 0.3563 

0.3559 0.3561 

0.3558 0.3560 

0.3557 0.3559 

0.3556 0.3557 

0.3554 0.3556 

0.3553 0.3555 



$iv(l) AN-'^ 

2.3779 2.3779 

2.3780 2.3780 

2.3781 2.3781 

2.3782 2.3782 

2.3783 2.3783 

2.3784 2.3784 

2.3785 2.3785 

2.3786 2.3786 

2.3787 2.3787 

2.3788 2.3787 







AN--^ 


0. 


,9688 


0, 


.9689 


0. 


,9687 


0, 


.9687 


0. 


,9685 


0, 


.9686 


0. 


,9684 


0, 


.9684 


0. 


,9682 


0, 


.9683 


0. 


,9680 


0, 


.9681 


0. 


,9679 


0, 


.9679 


0. 


,9677 


0, 


.9678 


0. 


,9676 


0, 


.9676 


0. 


,9674 


0, 


.9675 


0. 


,9673 


0, 


.9673 


^ 


n{1) 


AN-'^ 


0. 


,3552 


0, 


.3553 


0. 


,3551 


0, 


.3552 


0. 


,3549 


0, 


.3551 


0. 


,3548 


0, 


.3550 


0. 


,3547 


0, 


.3548 


0. 


,3546 


0, 


.3547 


0. 


,3544 


0, 


.3546 


0. 


,3543 


0, 


.3544 


0. 


,3542 


0, 


.3543 


0. 


,3541 


0, 


.3542 


0. 


,3539 


0, 


.3540 



$^(1) AN- 



2.3789 


2, 


.3788 


2.3789 


2, 


.3789 


2.3790 


2, 


.3790 


2.3791 


2, 


.3791 


2.3792 


2, 


.3792 


2.3793 


2, 


.3793 


2.3794 


2, 


.3794 


2.3795 


2, 


.3795 


2.3796 


2, 


.3796 


2.3797 


2, 


.3797 


$iv(l) 


AN-"^ 


0.9671 


0, 


.9672 


0.9670 


0, 


.9670 


0.9668 


0, 


.9669 


0.9667 


0, 


.9667 


0.9665 


0, 


.9665 


0.9664 


0, 


.9664 


0.9662 


0, 


.9662 


0.9661 


0, 


.9661 


0.9659 


0, 


.9659 


0.9658 


0, 


.9658 


0.9656 


0, 


.9656 


$7V(1) 


AN-'^ 


0.3538 


0, 


.3539 


0.3537 


0, 


.3538 


0.3536 


0, 


.3537 


0.3535 


0, 


.3535 


0.3533 


0, 


.3534 


0.3532 


0, 


.3533 


0.3531 


0, 


.3532 


0.3530 


0, 


.3530 


0.3528 


0, 


.3529 


0.3527 


0, 


.3528 


0.3526 


0, 


.3526 



#jv(l) 


AN-'^ 


2.3798 


2, 


.3798 


2.3799 


2, 


,3799 


2.3800 


2, 


.3800 


2.3801 


2, 


.3801 


2.3802 


2, 


.3802 


2.3803 


2, 


,3803 


2.3804 


2, 


.3804 


2.3805 


2, 


.3805 


2.3805 


2, 


.3805 


2.3806 


2, 


.3806 



#iv(l) AN-'^ 

0.9654 0.9655 

0.9653 0.9653 

0.9651 0.9652 

0.9650 0.9650 

0.9648 0.9649 

0.9647 0.9647 

0.9646 0.9646 

0.9644 0.9644 



$iv(l) AN-<^ 

0.3525 0.3525 

0.3524 0.3524 

0.3522 0.3523 

0.3521 0.3521 

0.3520 0.3520 

0.3519 0.3519 

0.3518 0.3518 

0.3517 0.3517 
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